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METHOD AND APPARATUS FOR GENERATING 
TEMPORALLY INTERPOLATED PROJECTIONS 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional Application No. 
60/397,658 filed on July 23, 2002 and U.S. Provisional Application 60/398,463 filed 
on July 25,2002. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to the field of medical imaging and 
more specifically to the field of imaging dynamic, internal tissue, such as cardiac 
tissue, by computed tomography. In particular, the present invention relates to the 
generation and reconstruction of temporally interpolated projection data. 

Computed tomography (CT) imaging systems measure the attenuation of X-ray 
beams passed through a patient from numerous angles. Based upon these 
measurements, a computer is able to reconstruct images of the portions of a patient's 
body responsible for the radiation attenuation. As will be appreciated by those skilled in 
the art, these images are based upon separate examination of a series of angularly 
displaced images of the transmitted X-ray beam. A CT system processes X-ray intensity 
data to generate 2D maps of the line integral of linear attenuation coefficients of the 
scanned object at multiple view angle positions, denoted as projection data. These data 
are then reconstructed to produce an image, which is typically displayed on a monitor, 
and may be printed or reproduced on film. A virtual 3-D image may also be produced 
by a CT examination. 

CT scanners operate by projecting fan shaped or cone shaped X-ray beams 
from an X-ray source. The X-ray beams may be collimated to control the shape and 
spread of the beams. The X-ray beams are attenuated as they pass through the object 
to be imaged, such as a patient. The attenuated beams are detected by a set of detector 
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elements. Each detector element produces a signal affected by the attenuation of the 
X-ray beams, and the data are processed to produce signals that represent the line 
integrals of the attenuation coefficients of the object along the X-ray paths. These 
signals are typically called "projection data" or just "projections". By using 
5 reconstruction techniques, such as filtered backprojection, useful images may be 

formulated from the projections. The images may in turn be associated to form a 
volume rendering of a region of interest. The locations of regions of interest, such as 
pathologies, may then be identified either automatically, such as by a computer- 
assisted detection (CAD) algorithm or, more conventionally, such as by a trained 
10 radiologist. CT scanning provides certain advantages over other types of techniques in 

diagnosing disease, particularly because it illustrates the accurate anatomical 
information about the body. Further, CT scans may help physicians distinguish between 
types of abnormalities more accurately. 

15 CT imaging techniques, however, may present certain challenges when imaging 

dynamic internal tissues, such as the heart. For example, in cardiac imaging, the motion 
of the heart causes inconsistencies in the projection data, which, after reconstruction, 
may result in various motion-related image artifacts such as blurring, streaking, or 
discontinuities. To reduce the occurrence of motion-related image artifacts, various 

20 techniques may be employed to improve the temporal resolution of the imaging system, 

thereby reducing the effects of the moving tissue. Temporal resolution may generally be 
improved by decreasing the rotation time of the CT gantry. In this way, the amount of 
motion that occurs within the temporal window associated with the acquisition of a 
projection data set is minimized. 

25 

Temporal resolution may be further improved by the choice of reconstruction 
algorithm. For example, segment reconstruction algorithms, such as half-scan 
reconstruction algorithms, may be employed in the reconstruction process. The segment 
reconstruction algorithms typically reconstruct images using projection data collected 
30 over an angular range of 180° plus the fan angle (a) of the X-ray beam. Because the 

acquisition of projection data during gantry rotation of 180° + a requires less time when 
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compared to acquisition during 360° of gantry rotation, the temporal resolution in the 
reconstructed images is improved. 

Multi-sector reconstruction techniques may also improve the temporal resolution 
5 of the reconstructed images by using projection data acquired during multiple rotations 

of the gantry by a multi-slice detector array. The projection data set used for 
reconstruction is composed of two or more sectors of projection data that are acquired 
from different cardiac cycles. The sectors comprise the data acquired during a short 
span of the gantry rotation, typically less than half of a rotation. The sectors, therefore, 
1 0 have good temporal resolution if acquired by a rapidly rotating gantry, thereby providing 

a good effective temporal resolution for the aggregate projection data set used in 
reconstruction. 

Using the techniques discussed above, third and fourth generation CT systems 
15 are capable of temporal resolutions of approximately 300 ms using segment 

reconstruction techniques. Fifth generation CT systems, utilizing a stationary detector 
ring and an electron gun which sweeps an electron beam across a stationary target ring 
to generate x-rays, are capable of achieving a temporal resolution of approximately 50 
ms. A temporal resolution of approximately 20 ms, however, is desirable in order to 
20 "freeze" cardiac motion, thereby minimizing motion related artifacts in the reconstructed 

images. For third generation CT systems, improving temporal resolution by the above 
techniques has typically focused on further increasing the rotational speed of the gantry. 

However, as the rotational speed of the gantry increases, the centripetal force on 
25 the gantry components also increases. The increasing centripetal force and the 

tolerances of the gantry components may comprise, therefore, a mechanical limitation to 
increases in gantry angular velocity. Furthermore, to obtain consistent image quality in 
terms of signal-to-noise ratio, a constant X-ray flux should be delivered to the imaged 
object or patient during the scan interval. Achieving a constant X-ray flux, however, 
30 places increased demand on the X-ray tube, particularly in regard to tube output, and on 

the components that cool the X-ray tube. Both mechanical and X-ray flux 
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considerations, therefore, are obstacles to increasing the gantry rotation speed 
sufficiently to achieve a temporal resolution of 20 ms or better in CT reconstructions. A 
technique for achieving a temporal resolution without increasing gantry rotation speed is 
therefore desirable. 

5 

BRIEF DESCRIPTION OF THE INVENTION 

The present technique provides a novel method and apparatus for improving 
temporal resolution of a CT imaging system. The technique employs a slowly rotating 
distributed X-ray source used to acquire projection data at each view position of an 
10 object or patient. The projection data, based upon the cardiac phase data and the 

frequency content of the projection data, may be filtered and temporally interpolated 
to generate interpolated projection data corresponding to a particular instant in the 
cardiac cycle. A high temporal resolution CT image may be reconstructed from the 
interpolated projection data. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other advantages and features of the invention will become 
apparent upon reading the following detailed description and upon reference to the 
drawings in which: 

20 

Fig. 1 is a diagrammatical view of an exemplary imaging system in the form of 
a CT imaging system for use in producing processed images, in accordance with one 
aspect of the present technique; 

25 Fig. 2 is another diagrammatical view of a physical implementation of the CT 

system of Fig. 1, in accordance with one aspect of the present technique; 

Fig. 3 is a diagrammatical view of an exemplary CT gantry, in accordance with 
one aspect of the present technique; 

30 
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Fig. 4 is a flowchart depicting the acquisition of a projection data set, in 
accordance with one aspect of the present technique; 

Fig. 5 is a sequence of views of the operation of a distributed X-ray source 
5 over the course of two cardiac cycles, in accordance with one aspect of the present 

technique; and 

Fig. 6 is a flowchart depicting the generation and reconstruction of an 
interpolated projection data set. 

10 

DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS 

Fig. 1 illustrates diagrammatically an imaging system 10 for acquiring and 
processing image data. In the illustrated embodiment, system 10 is a computed 
tomography (CT) system designed to acquire X-ray projection data, to reconstruct the 
15 projection data into an image, and to process the image data for display and analysis in 

accordance with the present technique. In the embodiment illustrated in Fig. 1, 
imaging system 10 includes a distributed source of X-ray radiation 12 positioned 
adjacent to a collimator 14. In this exemplary embodiment, the distributed source of 
X-ray radiation 12 comprises one or more addressable X-ray focal spots. 

20 

Collimator 14 permits a stream of radiation 16 to pass into a region in which a 
subject, such as a human patient 18 is positioned. The stream of radiation 16 may be 
generally cone shaped, depending on the configuration of the detector array, discussed 
below, as well as the desired method of data acquisition. A portion of the radiation 20 
25 passes through or around the subject and impacts a detector array, represented 

generally at reference numeral 22. Detector elements of the array produce electrical 
signals that represent the intensity of the incident X-ray beam. These signals are 
acquired and processed to reconstruct an image of the features within the subject. 

30 The distributed X-ray source 12 is controlled by a system controller 24, which 

furnishes both power, focal spot location, and control signals for CT examination 
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sequences. Moreover, detector 22 is coupled to the system controller 24, which 
commands acquisition of the signals generated in the detector 22. The system 
controller 24 may also execute various signal processing and filtration functions, such 
as for initial adjustment of dynamic ranges, interleaving of digital image data, and so 
5 forth. In general, system controller 24 commands operation of the imaging system to 

execute examination protocols and to process acquired data. In the present context, 
system controller 24 also includes signal processing circuitry, typically based upon a 
general purpose or application-specific digital computer, associated memory circuitry 
for storing programs and routines executed by the computer, as well as configuration 
10 parameters and image data, interface circuits, and so forth. 

In the embodiment illustrated in Fig. 1, system controller 24 is coupled to a 
linear positioning subsystem 26 and rotational subsystem 28. The rotational 
subsystem 28 enables the distributed X-ray source 12, collimator 14 and the detector 

15 22 to be rotated one or multiple turns around the patient 18. It should be noted that 

the rotational subsystem 28 might include a gantry. Thus, the system controller 24 
may be utilized to operate the gantry. The linear positioning subsystem 26 enables the 
patient 18, or more specifically a patient table, to be displaced linearly. Thus, the 
patient table may be linearly moved within the gantry to generate images of particular 

20 areas of the patient 18. Although the embodiment of the system described above is a 

modified third generation CT scanner, in which the distributed source 12 and the 
detector 22 both rotate, the methods to generate signals representative of cardiac 
motion described herein apply to all advanced generation CT systems. 

25 As will be appreciated by those skilled in the art, the distributed source of 

radiation may be controlled by an X-ray controller 30 disposed within the system 
controller 24. Particularly, the X-ray controller 30 is configured to provide power, 
focal spot location, and timing signals to the distributed X-ray source 12. A motor 
controller 32 may be utilized to control the movement of the rotational subsystem 28 

30 and the linear positioning subsystem 26. 
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Further, the system controller 24 is also illustrated comprising a data 
acquisition system 34. In this exemplary embodiment, the detector 22 is coupled to 
the system controller 24, and more particularly to the data acquisition system 34. The 
data acquisition system 34 receives data collected by readout electronics of the 
5 detector 22. The data acquisition system 34 typically receives sampled analog signals 

from the detector 22 and converts the data to digital signals for subsequent processing 
by a computer 36. 

The computer 36 is typically coupled to the system controller 24. The data 
10 collected by the data acquisition system 34 may be transmitted to the computer 36 for 

subsequent processing and reconstruction. The computer 36 may comprise or 
communicate with a memory 38 that can store data processed by the computer 36 or 
data to be processed by the computer 36. It should be understood that any type of 
computer accessible memory device capable of storing the desired amount of data 
15 and/or code may be utilized by such an exemplary system 10. Moreover, the memory 

38 may comprise one or more memory devices, such as magnetic or optical devices, 
of similar or different types, which may be local and/or remote to the system 10. The 
memory 38 may store data, processing parameters, and/or computer programs 
comprising one or more routines for performing the processes described herein. 

20 

The computer 36 may also be adapted to control features enabled by the 
system controller 24, i.e., scanning operations and data acquisition. Furthermore, the 
computer 36 may be configured to receive commands and scanning parameters from 
an operator via an operator workstation 40 typically equipped with a keyboard and 
25 other input devices (not shown). An operator may thereby control the system 10 via 

the input devices. Thus, the operator may observe the reconstructed image and other 
data relevant to the system from computer 36, initiate imaging, and so forth. 

A display 42 coupled to the operator workstation 40 may be utilized to observe 
30 the reconstructed image. Additionally, the scanned image may also be printed by a 

printer 44 which may be coupled to the operator workstation 40. The display 42 and 
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printer 44 may also be connected to the computer 36, either directly or via the operator 
workstation 40. Further, the operator workstation 40 may also be coupled to a picture 
archiving and communications system (PACS) 46. It should be noted that PACS 46 
might be coupled to a remote system 48, radiology department information system 
(RIS), hospital information system (HIS) or to an internal or external network, so that 
others at different locations may gain access to the image data. 

It should be further noted that the computer 36 and operator workstation 40 
may be coupled to other output devices, which may include standard, or special 
purpose computer monitors and associated processing circuitry. One or more operator 
workstations 40 may be further linked in the system for outputting system parameters, 
requesting examinations, viewing images, and so forth. In general, displays, printers, 
workstations, and similar devices supplied within the system may be local to the data 
acquisition components, or may be remote from these components, such as elsewhere 
within an institution or hospital, or in an entirely different location, linked to the 
image acquisition system via one or more configurable networks, such as the Internet, 
virtual private networks, and so forth. 

Referring generally to Fig. 2, an exemplary imaging system utilized in a 
present embodiment may be a CT scanning system 50. The CT scanning system 50 
may be a volumetric CT (VCT) system utilizing a cone-beam geometry and an area 
detector to allow the imaging of a volume, such as an entire internal organ. 
Furthermore, as noted above, the CT scanning system 50 may be a modified third 
generation CT imaging system, as depicted, or may be a later generation CT imaging 
system. 

The CT scanning system 50 is illustrated with a frame 52 and a gantry 54 that 
has an aperture 56 through which a patient 18 may be moved. A patient table 58 may 
be positioned in the aperture 56 of the frame 52 and the gantry 54 to facilitate 
movement of the patient 18, such as via linear displacement of the table 58 by the 
linear positioning subsystem 26 (see Fig. 1). The gantry 54 is illustrated with a 
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distributed source of radiation 12, typically one or more X-ray sources that emit X-ray 
radiation from one or more focal points 62. The distributed X-ray source 12 may 
change the location of a focal point 62 instantaneously, either longitudinally or along 
the arc of the gantry 54. For cardiac imaging, the stream of radiation 16 is directed 
5 towards a cross section of the patient 18 including the heart. 

In typical operation, distributed X-ray source 12 projects X-rays from the one 
or more focal points 62 and toward the detector array 22. The collimator 14 (see Fig. 
1), such as lead or tungsten shutters, typically defines the size and shape of the 

10 projected X-rays that emerge from the distributed X-ray source 12. The detector 22, 

such as an area detector in the case of a VCT system, is generally formed by a 
plurality of detector elements, which detect the X-rays that pass through and around a 
subject of interest, such as the heart or chest. Each detector element produces an 
electrical signal that represents the intensity of the X-ray beam at the position of the 

1 5 element during the time the beam strikes the detector. The gantry 54 is rotated around 

the subject of interest so that a plurality of radiographic views may be collected by the 
computer 36. 

Thus, as the distributed X-ray source 12 and the detector 22 rotate, the detector 
20 22 collects data of the attenuated X-ray beams. Data collected from the detector 22 

then undergoes pre-processing and calibration to condition the data to represent the 
line integrals of the attenuation coefficients of the scanned objects. The processed 
data, commonly called projections, may then be filtered and backprojected to 
formulate an image of the scanned area. A formulated image may incorporate, in 
25 certain modes, projection data acquired for an angular rotation of the gantry of less or 

more than 360 degrees. 

Once reconstructed, the cardiac image produced by the system of Figs. 1 and 2 
reveals the heart of the patient 18. As illustrated generally in Fig. 2, the image 64 may 
30 be displayed to show patient features, such as indicated at reference numeral 66 in Fig. 

2. In traditional approaches to diagnosis of medical conditions, such as disease states, 
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and more generally of medical conditions or events, a radiologist or physician would 
consider the reconstructed image 64 to discern characteristic features of interest. Such 
features 66 include coronary arteries or stenotic lesions of interest, and other features, 
which would be discernable in the image, based upon the skill and knowledge of the 
5 individual practitioner. Other analyses may be based upon capabilities of various 

CAD algorithms. 

As will be appreciated by those skilled in the art, reconstruction of an image 
64 may be complicated by a variety of factors. For example, reconstructed images 64 
10 of dynamic tissue may include motion-related image artifacts that are attributable to 

the movement of the tissue during imaging. To reduce motion-related artifacts, it is 
generally desirable to improve the temporal resolution of the CT reconstruction 
process. 

To address this problem, the present technique uses a projection data set, along 
with frequency content and phase information associated with the projection data set, 
to generate a set of temporally interpolated projections. The temporally interpolated 
projections may be reconstructed to form images with a high temporal resolution, 
typically less than 50 ms. One approach of the present technique utilizes a CT gantry 
54 in which the scan path of the gantry 54 is divided into a discrete number, N, of 
uniformly sized arcs 80, as depicted in Fig. 3. For example, as depicted in Fig. 3, N 
equals 18. The distributed X-ray source 12 may be sized to cover two arc lengths 82. 
Alternatively, the distributed source 12 may cover less than or more than two arc 
lengths if alternative sampling schemes are employed during acquisition of the 
radiographs. 

In the present technique, as depicted at step 90 of Fig. 4, the distributed X-ray 
source 12 is slowly rotated about a subject to generate radiographs, which may 
subsequently be processed into interpolated projection data. For example, the 
30 distributed X-ray source 12 may execute a single rotation about the subject, during 

which multiple cardiac cycles may occur. The rotation rate of the distributed source 
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X-ray source 12, and the detector 22 in suitably modified third generation CT systems, 
may be greater than eight seconds, and may be between ten seconds and thirty seconds 
in one implementation. In particular, the rotation rate may equal the time required for 
one cardiac cycle of the subject multiplied by the number of arc lengths, N. 

5 

Projection data is repeatedly acquired as the distributed source 12 passes along 
each arc 80, as depicted at step 92. If the rotation rate of the distributed source 12 
equals the time required for one cardiac cycle multiplied by N 9 the start of a new 
cardiac cycle should coincide with the start of projection acquisition along the next arc 
10 80. The position of the focal point 62 on the distributed X-ray source 12 may be 

continually adjusted as the gantry 54 rotates so that each sampled view position 
relative to the object being scanned is identical. 

For example, referring to Fig. 5, the rotation and acquisition process at two 
15 sequential arcs 80 is depicted. To simplify presentation, the cardiac cycle is presumed 

to be one second. At the outset, i.e., T = 0 seconds, the leading edge 100 of the 
distributed X-ray source 12 coincides with the boundary between the first arc 102 and 
the second arc 104, such that the first arc 102 is covered by half of the distributed X- 
ray source 12. The portion of the distributed X-ray source 12 coinciding with the first 
20 arc 102 emits X-rays 16 through the volume of the thorax including the subject's heart 

106 by triggering the focal spots along the arc in some predefined order. The portion 
of the distributed X-ray source 12 that has yet to enter the first arc 102 does not emit 
X-rays 16. 

25 At T = 0.5 seconds, the distributed X-ray source 12 has rotated halfway 

through the first arc 102 such that the middle portion of the distributed X-ray source 
12 coincides with the first arc 102. The portion of the distributed X-ray source 12 
coinciding with the first arc 102 continues to emit X-rays 16 through the volume of 
the thorax occupied by the subject's heart 106. That portion of the distributed X-ray 

30 source 12 which has not yet entered the first arc 102 or which has passed into the 

second arc 104 does not emit X-rays 16. 
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At T= 1.0" seconds, i.e., immediately before 1.0 seconds and the start of the 
next cardiac cycle, the trailing edge 108 of the distributed X-ray source 12 coincides 
with the boundary between the first arc 102 and the preceding arc 1 10. The portion of 
the distributed X-ray source 12 coinciding with the first arc 102 continues to emit X- 
rays 16 through the volume of the thorax occupied by the subject's heart 106. That 
portion of the distributed X-ray source 12 which has passed into the second arc 104 
does not emit X-rays 16. 

Because the X-ray source 12 is distributed, the switching of the arc positions at 
T = 1.0 seconds occurs instantaneously. As a result, at T = 1.0 + seconds, i.e., 
immediately after 1.0 seconds, the leading edge 100 of the distributed X-ray source 12 
coincides with the boundary between the second arc 104 and the subsequent arc 112. 
The acquisition of projections along the second arc 104 thereafter precedes as 
described with regard to the first arc 102, as depicted at T= 1.5 seconds and 7= 2.0" 
seconds. In this manner, projection data is acquired along each arc 80 repeatedly 
during the respective cardiac cycle, as depicted at step 92 in Fig. 4. If the rotation or 
rotations are not complete, projection data may be acquired for the next arc 80 and 
cardiac cycle, as incremented at step 116, until the rotation or rotations are complete, 
as determined at decision block 118. The result is a projection data set 120 that 
comprises projections acquired at different instants in time with respect to the cardiac 
cycle at each view position. 

Although the distributed X-ray source 12 is shown to cover two arc segments 
80 in Fig. 3 and in Fig. 5, the distributed source 12 can be made slightly larger than 
this to account for variations in the cardiac cycle during the data acquisition process. 
Alternatively, the distributed X-ray source 12 may be sized to cover one arc segment 
80, or slightly larger than one arc segment 80 to account for variations in the cardiac 
cycle during data acquisition, and used in conjunction with alternative sampling 
routines to achieve the same effective sampling data set. 
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The frequency content 130 of the projection data, i.e., a priori information 
about characteristics of motion in the heart as represented in the signals measured by 
elements of detector 22, and the cardiac phase data 132, i.e., the timing of the cardiac 
phases during the acquisition of the projection data 120, may be used to interpolate the 
5 projection data 120, as depicted in Fig. 6. The phase data 132 may comprise an ECG 

signal acquired concurrently with the radiographs. Alternatively, phase data 132 may 
be derived from the projection data 120, such as via techniques employing consistency 
conditions to analyze the projection data 120 and/or to compare the moments of the 
projection data 120. 

10 

Using the phase data 132 and information about the frequency content 130 in 
the projection data, the projections are correlated with the times that they correspond 
to in the cardiac cycle as well as with the angular positions of the gantry 54 to which 
they correspond. The projection data 120 may then be interpolated, as depicted at step 

15 134, to generate interpolated projections 136. Because the two-dimensional 

projections appropriately capture the frequency content 130 in the signal acquired at 
each detector element, the interpolated projections 136 each correctly characterize the 
projection data at any instant in time with respect to the cardiac cycle at the respective 
view location. In this manner, projections acquired at discrete points in time may be 

20 converted into a continuous-time representation and, from the continuous time 

representation, the projections may be interpolated to a particular instant with respect 
to the cardiac cycle. 

The conversion of discrete points into a continuous-time representation may be 
25 accomplished in a variety of ways known in the art. Likewise, the interpolation of 

values from a continuous-time representation using a suitable interpolation algorithm 
may be accomplished in various ways. For example, the Nyquist Theorem provides 
that, if a sufficiently high rate of discrete samples of a waveform are obtained, a 
continuous-time representation using the discrete-time samples may be generated. 
30 The Nyquist Theorem also provides that the value of any sample along the 

continuous-time representation can be interpolated at any point in time by selecting 
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the sample value on the continuous-time representation that corresponds to the point 
in time. For example, a Fourier time series is a suitable continuous-time function for 
this purpose if the motion is periodic. 

As noted above, the interpolated projections 136 each correspond to a 
particular instant of the cardiac cycle at a respective view angle. The interpolated 
projections 136 corresponding to a desired instant of the cardiac cycle may therefore 
be reconstructed, as depicted at step 138, to generate cardiac images at the desired 
instant of the cardiac cycle. The reconstructed images may, if desired be associated 
spatially and/or temporally to generate an image over time, a volume 140 at an instant 
in time, or a volume 140 over time. In addition, because the interpolated projections 
136 are interpolated to the same instant in time, the reconstructed images and/or 
volumes have a high temporal resolution, typically less than 50 ms. Therefore the 
reconstructed images and/or volumes are substantially free of motion defects and 
artifacts, effectively "freezing" the cardiac motion at each point in time. 

In addition, the interpolation process provides a mechanism to reduce quantum 
noise in the projection data. For example, a priori information about the frequency 
content of the projection data 120 may be used to bandlimit frequency components in 
20 the projection data 120. Bandlimiting the frequency components may help reduce 

noise in the reconstructed images and may allow patient dose to be reduced while still 
achieving suitable quality of reconstructed images. 

While the present techniques have been presented in the context of cardiac 
25 imaging, the techniques may be applied to the imaging of other dynamic objects. 

Discussion of cardiac imaging is presented merely to facilitate explanation of the 
present techniques. Indeed, while the invention may be susceptible to various 
modifications and alternative forms, specific embodiments have been shown by way 
of example in the drawings and have been described in detail herein. However, it 
30 should be understood that the invention is not intended to be limited to the particular 

forms disclosed. Rather, the invention is to cover all modifications, equivalents, and 
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alternatives falling within the spirit and scope of the invention as defined 
following appended claims. 
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